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Abstract 

Using Quantum Field Theory we derive a general formula for the double inclusive 
spectra of photons radiated by a system in local equilibrium. The derived expres- 
sion differs significantly from the one mostly used up to now in photon intensity 
interferometry of heavy-ion collisions. We present a covariant expression for dou- 
ble inclusive spectra adapted for usage in numerical simulations. Application to 
a schematic model with a Bj0rken type expansion gives strong evidence for the 
need of reinvestigating photon-photon correlations for expanding sources. 
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Double inclusive spectra of hard photons radiated at the early stage of heavy-ion col- 
lisions have become an important subject of theoretical [1, 2, 3] and experimental [4] 
investigations. Since these photons originate from very energetic collisions and interact 
very little with the matter they are traversing, they may carry signatures of the hot and 
dense phase of the matter, in particular of the quark-gluon plasma. 

Except for the work of D. Neuhauser [1] in which he treated a static source, usually one 
is analyzing the correlations in two-photon coincidence measurements with expressions 
based on the corresponding ones for scalar bosons [5] augmented with a degeneracy factor 
which is supposed to take care of the spin structure of the photons. 

The aim of this work is to derive in a quantum-field treatment covariant expressions for 
the two-photon density and herewith the corresponding correlation function. In particular 
we pay attention to charge conservation and to the fact that massless photons have only 
two helicity eigenstates. Due to these properties averaging over all spin directions results 
in a reduction factor of one-half compared to the scalar case only if the momenta of 
the two photons are parallel. For a finite angle, however, the interference is reduced, so 
that the deduced size of the source is smaller than the one obtained from the correlation 
function for scalar bosons. For Bj0rken type hydrodynamics we estimate a difference in 
the apparent source size of about 30% in rapidity. Therefore, we conclude that the "pion- 
inspired" formulas for the Bose-Einstein correlation function is not applicable for photon 
intensity interferometry studies, especially if one deals with an expanding source. 

The finally measured number of photons per invariant momentum bin with polariza- 
tion A is given by the density matrix of the final state pf which can be expressed in terms 
of the initial density matrix pi and the S'-matrix using the relation Pf — SpiS^. 

PlX{k) = k °^k = MP/4(fc)cA(fc)} = Tr{p^c[(k)c x (k)S} (1) 
The same consideration as above leads to the expression for the number of photon pairs 

M l\jXl\2 

Pt X2 {k,M) = k i k 2 d 3 kid s k2 = rr{p / 4 2 (fc 2 )4 1 (fc 1 )c Al (fc 1 )c A2 (fc 2 )} = 

= Tr{p^t4 2 (fc 2 )4 i (fc 1 )cA 1 (fc 1 )c A2 (fc 2 )5} . (2) 

If we separate the interaction part of the Lagrangian into for strong interactions 

and J ft (x)A li (x) for the coupling to the electromagnetic field the S'-matix can be written 
in the interaction picture by the time-ordered exponential [6] 

S = T exp J d 4 x(Li nt (x) + J^{x)A tl (x))^ . (3) 
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To proceed further we assume [7] that the initial state does not contain hard photons, 
which means 

cx{k) Pl = Pi c\{k) = . (4) 
The condition (4) together with the identities 

'c x (k),S]=iT(j x (k)S) , (5) 

'_S\c[(k)] =-iT(S*j1(kj) , 

where = (\k\,k) and J\(k) is the Fourier transform of the transverse current operator 
given by 

J x {k)^ei{k)Jd'xe^ k r{x) (6) 

J x \k)^ei(k)Jd 4 xe-^r(x) , 

allow us to write one- and two-body densities via the chronological (T) and antichrono- 
logical (T) products of the current operators: 

P*(fc) = Tr{ Pi f(tfJ x \k))T(j x (k)S)} (7) 

P 2 AlAa (fei,fc 2 ) =Tr| Pt t(5 t J A t 2 (A; 2 )J A t i (A; 1 ))T(j Al (A; 1 )JA 2 (A: 2 )5)} . (8) 

All currents in the expressions written above are operators in the interaction representa- 
tion. Utilizing properties of time ordering and of the S'-matrix we can write the one- and 
two-photon spectra in terms of the current operators in the Heisenberg representation 
(J) [8] as 

P?(k) =Tri yPl j x \k)Uk)) (9) 

P 2 AlAa (fci,fc 2 ) =Tr{p l t(j A t 2 (A; 2 )J A t i (A; 1 ))T(j Al (A; 1 )JA 2 (A: 2 ))} . (10) 

Up to now we have made no approximation to tackle the unsolvable many-body problem. 
Expressions (9) and (10) are still exact but they are recast into a form which is better 
suited for approximations than the original definitions (1) and (2). Formally they contain 
only operators of the strongly interacting system. 

In the following we apply eqs. (9) and (10) to a highly excited system which due to 
short-ranged hard collisions of the charged constituents is radiating energetic photons. 
The collisions cause locally rapid changes in the electric charge current density J^(x) , 
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which result in non-zero Fourier components in J^(k) for large \k\ . The energetic photons 
are weakly interacting so that, once they are created, they are assumed to leave the 
hadron system without further interactions. Therefore, for calculating the final photon 
distributions one needs a local production rate which is then integrated over space and 
time. Eq. (9) is of this type. In coordinate space it reads 

Pi\k) = e^(fc)eJ 2 (A;)/^ 4 Axe-^rr{p,>(x + ^)> 2 (x - ^)} (11) 

where the local rate is given by the integral over Ax of the current-current correlator. 

The following approximations, which will lead to an expression for the correlator 
adapted for usage in numerical simulations for heavy-ion collisions, are based upon three 
physical assumptions: 

(i) The hadronic correlations in the system are of short range in space-time. The typical 
correlation length is about the mean free path \ s of the strongly interacting particles. 
Therefore, 

<A(x)B(y)>xi<A(x)><B(y)> if \x - y\ > X s or \x°-y°\>\ s . 

(ii) The hadronic mean free path X s is much less than the characteristic size of variations 
of the macroscopic variables in space and time, denoted by L . 

(Hi) Only hard photons with high momenta \k\ are considered ( jfcj A s > 1 ) , for which the 
slowly varying collective current <J^(x)> (proportional to the collective 4- velocity 
field ^(x) ) does not contribute, i.e. <J^(k)>^ for high enough \k\ . 

For instance, applying these assumptions to (11) means that the integral in Ax effectively 

spreads over the 4-volume of the order Af and the single inclusive spectrum -Pi(fc) is 

therefore proportional to L 4 X^ . 

Inside the trace on the right-hand side of (11) the mean time x° can be moved over 

to the statistical operator such that 

Tr{pJ^(x + ^)>(x - ^)} = Tr{p(x°)J^(x + ^)>(5 - ^)} , (12) 

where x = (0, x) and p(x°) = e~ lHx ° pie lHx ° is the solution of the Liouville equation for 
the density matrix. For the case of local equilibrium, where the time-dependence enters 
the density matrix p(x°) only through the thermodynamic quantities, like for example 
temperature T(x°,x) , 4-velocity u^(x°,x) etc., the identity (12) is a starting point to 
express the current-current correlator in terms of thermodynamic variables. 
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The double inclusive spectrum is given by the four-point correlator (10) and like any 
other Green function can be written in terms of connected parts as 

P 2 AlAa (fci,fc 2 ) =P Al (fc 1 )P A2 (fc 2 ) + |<J1 1 (A; 1 )J A2 (A; 2 )>| 2 + \<T(j Xl (h) Jx 2 (h))>\ 2 

+ «t(j A t 2 (A; 2 )J A t i (A; 1 ))T(j Al (A; 1 )J A2 (A; 2 ))» . (13) 

All terms in (13) involve four space-time integrations. In coordinate space a connected 
correlator disappears if the distance between any pair of its space-time points exceeds 
the correlation length. The connected four-point correlator -C- O^ (last term on the 
right-hand side of (13) ) is only non-zero if the relative distance of all its space-time 
points are within X s . Therefore, this term is proportional to L 4 (A 4 ) 3 and hence, due to 
condition (ii) , is smaller by factor of (X s /L) 4 when compared to all other terms which are 
proportional to (L 4 ) 2 (A 4 ) 2 . The third term on the r.h.s of (13), \<T^J\ 1 (ki) </a 2 (^2))>| 2 , 
contains an integral Jd^x exp(i(|fci| + |fc 2 |)x°) • • • in which the energies of the two photons 
add up in the phase factor. This implies that only momenta which fulfill the condition 
(\ki\ + |fc 2 |) ■ L < 1 contribute appreciably to this term. These momenta are much smaller 
than required by conditions (ii) and (Hi) and hence we shall drop this term. 

After these approximations based on the statistical properties of the emitting source 
the one- and two-photon spectra (9) and (10) are given by the current-current correlator 
< J^^ki) J fJl2 (k 2 )> only. Because the electrical current is conserved (d lx J lx (x) = 0) the 
correlator must be transverse "from both sides" : 

fc Ml <> lt (A; 1 )> 2 (A; 2 )>=<> 1 t( A;i )> 2 ( A;2 ) > k ^ = . (14) 
In analogy to (11) the correlator can be written as 

<> lt (A; 1 )> 2 (A; 2 )>= ^e^^-^^^e^^-^^^-^^Axe^^ 1 ^ 2 )/ 2 

n 

xT r{ Pi J^(x + ^)J^(x-^)} , (15) 

where the whole space-time region occupied by the source is divided into cells with 4- 
volume Q n located at mean-coordinates x n . Under the assumptions (i) — (Hi) the size 
of the space-time cells should be chosen to be about \ s so that the radiation of hard 
photons from different cells can be considered as entirely independent processes. This 
implies that the currents are conserved for each cell separately and that the contribution 
of the cells to the total current correlator is additive. Therefore, we parametrize (15) as 

<> lt (fc 1 )>(fc 2 )> Y,^- l ^- k - i) Q^\k 1 M\Xn)w(k l M^n) , (16) 

n 
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where the tensor Q^ 1 ^ 2 carries the tensorial structure of the current-current correlator 
and the function w describes the strength. The Q-tensor is explicitly transverse 



h^Q^ik^hlx) = Q^ 2 (k u k 2 \x)k 2ll2 = (17) 

and, therefore, ensures current conservations for every cell. The tensor Q^ 1 ^ 2 and the func- 
tion w are hermitian in the sense that Q ,J,ltJ,2 (ki, k 2 \x)* = Q^ 12 ^ 1 (A; 2 , ki\x) and w(k±, k 2 \x)* = 
w(k 2 , ki\x). It is convenient to normalize the Q-tensor by Q£(k,k\x) = —2. From now 
on we will write an integral instead of the sum in (16). 

Altogether single and double inclusive photon spectra (no polarization measured) can 
be expressed with the help of (16) as follows 1 : 

2 r 

Px(fc) = Yl p i(k) = 2 / d 4 xw(k, k\x) (18) 

(A=l) 

P 2 (k 1 ,k 2 )= f; P 2 x ^ 2 (k 1 ,k 2 )=P 1 (k 1 )P 1 (k 2 ) (19) 

(Ai,A 2 =l) 

+ J Me-' M(tl " fe) »(fci, k 2 \x)w{k 2 , k^Q^ih, k 2 \x)Q^(k 2 , h\y) 

To find an explicit form of the Q-tensor we assume that the radiating system is in 
local equilibrium, that is, the radiation from an elementary cell is isotropic in its local 
rest frame and therefore Q ll%2 = S ll%2 (where i\,i 2 = 1,2,3). The other components 
of the Q-tensor are then uniquely determined by the transversality condition (17) as 
Q 00 = k ± k 2 /k^, Q il0 = k 2 x fk%, Q 0i2 = k i2 jk\. In local thermal equilibrium the only 
other 4- vector entering is the 4- velocity u^(x) of the emitting cell which is located at point 
x . Therefore, the only hermitian and transversal Q-tensor which can be constructed from 
ki , k 2 and u^(x) , and which in the rest frame of a cell (w M = (1, 0) ) reproduces the 
isotropic form given above, reads as follows: 

In eq. (20) we do not include the term proportional to k^k^ 2 because it does not con- 
tribute to observables and can be removed by an appropriate gauge transformation. 
1 We take the sum over photon polarisations by means of 

2 

X] e M ( fc ) e ^( fc ) = ~9iu> - kf,K/(sk) 2 + (k^s v + s^k v )/(sk) 

(A=l) 

( s M is some reference 4- vector s 2 = 1 ), where due to the transversality condition (17) only the first term 
on the right-hand side {—g^v) really contributes to (18), (19). 
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After defining the Q-tensor all dynamic information on the photon production is 
contained in the function w(ki, k 2 \x). Inserting (20) into (19) one gets 



P 2 (k l , k 2 )=P 1 {k 1 )P 1 {k 2 )+j d A xd% cos(AxAk)R{k u k 2 \x, y)w(h, k 2 \x)w{k 2 , h\y) , (21) 
where we introduce the abbreviation 

(u(x)k 2 )(u(y)ki) 
(fcifc 2 ) 1 (uWujytfjhfa) 2 , x 



The i?-function reflects the fact that photons are massless spin-1 particles and in this 
point expression (21) differs significantly from the one obtained in [5] for pions. The 
latter has been used for photons of one single polarization in [2, 3, 9, 10]. Since the 
experimental measurements will involve an averaging over polarizations, a comparison 
with data requires the new result (21). Being very sensitive to the 4-velocity field u^ix) 
the i?-function can change drastically the interference of two photons for a relativistically 
expanding source. But even for a static source the polarization properties of photons are 
very important. To illustrate this let us consider two photons with momentum k\ and 
k 2 which are emitted from two cells having the same 4-velocity u(x) — u(y) — (1, 0) . In 
that case [1] the i?-function reduces to 

2 2 

i?(A; 1; A; 2 |x,|/) = l + (cos^) 2 = £ (e Al (A:i)e A2 (A; 2 )) , (23) 

(Ai,A 2 =l) 

where cos9 = feife2/(|fei||fe2|) • The physical reason is easily understood for this special 
case when one considers the last part of equation (23). The summation over polarizations 
(here we choose the Coulomb gauge: e x (k) = (0,e x (k)) and ke x (k) = 0) does not just 
lead to a factor of two as implicitly assumed in [2, 3, 9, 10]. As illustrated by Fig. 1 only 
one direction of the linear polarization can be chosen equal for both photons, whereas the 
other polarization directions differ by the angle 9 between the two momenta. Therefore, 
the polarization overlap involves (cos#) 2 and is less than 2 for 9 ^ (see (23)). In 
a realistic situation the radiation comes of course from many cells with different four- 
velocities and the full structure of the Q-tensor has to be employed. Anyhow, for a 
non-vanishing angle between the observed photon momenta the correlation is reduced 
compared to the "pion inspired" recipe (25). 

For application in dynamical models we recommend to use the following formula: 

P 2 {k u k 2 ) =Pi(fci)Pi(fc 2 )+<> lt (A;i)> 2 (A; 2 )><j; 2 (A; 2 )J m (A; 1 )> (24) 
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Fig. 1: The linear polarization vectors e x (k) for two photons with momenta fci and ■ 

and calculate the current-current correlator with (16) keeping the explicit form of the 
Q-tensor as given in (20). This is equivalent to (21) but involves only one dh: integration 
instead of d^xdSj in (21). 

The dependence of w(ki,k 2 \x) on Ak = k\ — k 2 and k — (k± + k 2 )/2 has to be 
derived from a microscopic model for the photon production. The only general statement 
which can be made at this point is that due to condition (i) the current-current correlator 
decays in the local restframe as function of Ak on the long scale A^ 1 • In this paper we do 
not consider a specific microscopic model for w(ki, k 2 \x) but rather investigate the role of 
the tensorial structure of the correlator on photon-photon interferometry. Therefore we 
write w(ki, k 2 \x) = w(k, k\x) + 0(Ak^Ak v ) . For simplicity we drop terms of second order 
in Ak . The strength of photon production w(k, k\x) can to some extend be estimated 
from the single inclusive cross section [11]. 

Based on the discussion above we propose to reexamine the photon intensity interfer- 
ometry in the spirit of [2, 3, 9, 10] but using (24) for the double inclusive cross section. 

Unfortunately until now the Bose-Einstein correlations of photons have been studied 
using the formula for double inclusive spectra where the photons are assumed to be scalar 
massless particles, but normalizing the correlation function to 3/2 instead of 2 [9, 2, 10]. 
That formula is just expression (21) with R(k±, k 2 \x,y) = 2 and reads 

pwrang^ fc ^ = p^fc^p^) +2 J d^xdSj COs{AxAk) W {k, k\x)w{k, k\y) . (25) 

In order to demonstrate the size of the effect which arises from the vector nature of 
the photon and the conservation of electric charge we have performed numerical studies 
of photon intensity interferometry in Bj0rken hydrodynamics (equation of state p — e/3) 
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using both our formula (24) (or equivalently (21)) and the wrong one (25). We assume that 
the photons are produced from the expanding source in local equilibrium and parameterize 
the photon production rate as w(ki, k 2 \x) — w(k, k\x) = N- (T(x)) 2 exp(— (ku(x))/T(x)) 
(for a more precise expression cf. [11]). The plot of the correlation function defined as 
C 2 (ki,k 2 ) = P 2 (ki, k 2 ) / Pi(ki)Pi(k 2 ) shows a significant difference in the results based 
on the different formulas (see Fig. 2). The correct expression gives less correlations for 
non-zero relative rapidities which means that the source size deduced from experimental 
data is overestimated by the wrong formula (25). 

In summary we should like to stress that the derivation of a basic equation for the 
double inclusive spectrum is significantly modified by the fact that photons are massless 
particles with spin 1 and that they are produced by a conserved electric current. Un- 
der the conditions (i) — (in) which are fulfilled for high-energy photons in relativistic 
heavy-ion collisions, the expression (24) is a suitable starting point for photon intensity 
interferometry studies. 

We would like to acknowledge the fruitful discussions with G. Bertsch and J. Knoll 
as well as the comments of D. Seibert. Especially we are grateful to R. Weiner and 
M. Pliimer for critical remarks. 
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Fig. 2: Bose-Einstein correlation plot as a function of rapidity difference Ay at fixed transverse 
momenta = = 100 MeV for Bj0rken hydrodynamics with Tj = 200 MeV, 
Tf = 140 MeV and initial proper-time n — 0.3 fm/c . 

Solid line corresponds to the photon-photon correlations (24) (or equivalently (21)). 
Dashed line represents the wrong "pion-inspired" expression (25). 
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